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Abstract. As a follow up to articles dealing firstly with a convective variational 
formulation in a Milne-Cartan framework for non-dissipative multi fluid models, and 
secondly with various ensuing stress energy conservation laws and generalised virial 
theorems, this work continues a series showing how analytical procedures developed 
in the context of General Relativity can be usefully adapted for implementation in a 
purely Newtonian framework where they provide physical insights that are not so easy 
to obtain by the traditional approach based on a 3+1 space time decomposition. The 
present article describes the 4-dimensionally covariant treatment of various dissipa- 
tive mechanisms, including viscosity in non-superfluid constituents, superfluid vortex 
drag, ordinary resistivity (mutual friction) between relatively moving non-superfluid 
constituents, and the transvective dissipation that occurs when matter is transformed 
from one constituent to another due to chemical disequilibrium such as may be pro- 
duced by meridional circulation in neutron stars. The corresponding non dissipative 
limit cases of vortex pinning, convection, and chemical equilibrium are also considered. 
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1. Introduction 



This article continues the development [1, 2] of a coherent fully covariant approach 
to the construction and application of Newtonian fluid models of the more general kind 
required in the context of neutron star phenomena in cases for which it is necessary 
to allow for independent motion of neutronic and protonic constituents. Generalis- 
ing the approach that was originally introduced [3] for the special case of Landau's 
two-constituent superfluid model, the preceding articles dealt with idealised perfectly 
conservative models, for which a strictly variational formulation is available. In a com- 
plementary treatment by Prix [4], it has been shown how this covariant variational 
formulation can be translated into terms of the familiar kind of 3+1 direct product 
structure of space with time that has traditionally be used in a non-relativistic ap- 
proach, and also how it can be extended to allow for dissipative effects - such as chem- 
ical reactions and mutual resistivity between relatively moving currents of thermal and 
other kinds - while retaining as much as possible of the convenient [5, 6] mathematical 
machinery provided by the variational approach. In the same spirit, but continuing 
within a fully covariant framework, the purpose of the present article is to extend our 
treatment to allow for a wider range of dissipative mechanisms of the kind [7] likely 
to be relevant in neutron stars, particularly those due to the presence of superfluid 
vortices. 

In the non-dissipative applications considered in the preceding articles [1, 2] the 
status of any entropy current 4-vector s M that may have been involved was effectively 
the same as that of the other relevant conserved currents with 4- vectors n£ = n x ii£ for 
corresponding number density n£ and unit flow vectors u£ designated by various values 
of the chemical index label x. However, in the dissipative applications to be considered 
here, the entropy density s and the (no longer conserved) entropy current 



will have a privileged role, characterising a corresponding local thermal rest frame spec- 
ified by a unit 4-vector u£ for which the (barred) zero value, x=0, of the chemical index 
will be reserved, i.e. we shall set n% = s, = s M . 

While the other (particle) currents may still either be conserved, in the sense of 
having V^n^ = for certain values of the chemical index x, or else may have divergence 
V M n^ of unrestricted, positive or negative, sign - corresponding to the possibilities of 
particle creation or destruction - the second law of thermodynamics stipulates that 
entropy should never be destroyed, which means that in an entirely self contained 
treatment we must always have 



i.e. V^ni^ can never be negative. In some contexts [8] it may however be convenient 
to work in terms of an open (i.e. not completely self contained) model in which, al- 
though not actually destroyed, entropy is nevertheless effectively lost from the system 
by some local heat removal mechanism - such as the URCA (neutrino - antineutrino 
pair creation) process in a neutron star core - in which case the relevant remaining 




(1) 



V^>0, 
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entropy current s M would not necessarily have to respect the restriction (2), but would 
be subject to the modified inequality (23) that is given below. 

In the preceeding article [2] it was shown how, in a system governed by a multifluid 
action variation principle, the relevant (kinetic, internal, and gravitational) Lagrangian 
action contributions give rise to corresponding variationally defined stress energy con- 
tributions that combine to give a total T % v which satisfies a Noether type identity of 
the form 

V,T^=E/-, (3) 

in which, for each value of the chemical index X, the 4-covector denotes the varia- 
tionally defined non - gravitational force density acting on the correponding constituent. 
The strictly conservative case considered in the preceeding work was characterised by 
dynamical equations given, according to the variation principle, just by the requirement 
that each of the separate force densities should vanish. The purpose of the present 
article is to extend the analysis to a more general category of dynamical equations, 
whereby the force densities are not required to vanish but are given by non-conservative 
contributions from dissipative mechanisms of three different kinds, namely as viscos- 
ity, resistance against relative motion, and transfusion between the various chemical 
constituents. 

Although the suspension of the variation principle leaves a considerable amount of 
lattitude in the way the various kinds of dissipative force may be specified, the admis- 
sible forms of force law are considerably restricted by the requirement of compatibility 
with the second law of thermodynamics as embodied, for a self contained system, in 
the inequality (2). As in the preceding articles [1, 2] our work will be guided by pre- 
vious experience [9, 10, 11] with analogous dissipative effects in a General Relativistic 
framework, which (contrary to what is commonly supposed) is actually simpler for 
many purposes, and particularly for the treatment of electromagnetic effects, which 
are not included (except as possible external background forces) in the present strictly 
Newtonian analysis. 

2. Viscous stress 

The first of the dissipative mechanisms that we need to consider - and the only one 
that will occur in a single constituent fluid model - is that of viscosity, whose effect will 
be interpretable in terms of a gross stress energy density tensor 

X 

in which the total T ^ provided by the previously considered action contributions [2] 
is supplemented by further by viscous stress contributions that are not obtained 
from the Lagrangian action but that are included to allow for deviations (from what 
would otherwise be a local thermal equilibrium state) due to space gradients of the 
corresponding flow vectors u^. In accordance with what is suggested by detailed micro- 
scopic analysis of dilute gas models [12] it will be assumed that each such contribution 
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has contravariant components 

that are symmetric and purely spacelike, i.e. 

r x ^ = r x ^, r x "X = 0, (6) 

(where, as discussed in the preceeding work [1], 7 MP is the degenerate Newtonian space 
metric while t M is the preferred Newtonian time gradient) and it will also be assumed 
that the mixed version is strictly spacelike with respect to the corresponding fluid rest 
frame, i.e. for each value of the chemical index x we shall have 

<r^ = 0. (7) 

It follows that it will be expressible in the form 

^ M = 7 X ^ XPM , (8) 

where j xup is the positive indefinite (rank 3) space metric tensor that would be de- 
termined (in the manner described in the preceeding work [1]) by choosing the ether 
reference vector e M to coincide with the local flow vector u v , i.e. it is given by the 
defining relations 

7x^ = 7^, 7 X ,P< = 0, (9) 

with 

7& = (10) 

In order to set up an appropriate category of models, we proceed on the basis of 
the postulate that this gross stress energy tensor should satisfy an energy momentum 
balance condition of the form 

y rp fi Y fx (11) 

in which the terms on the right will all vanish whenever we are dealing with a strictly 
self contained system, but in which the possibility of external force density contributions 
/ e x is included to allow for cases when we are dealing with an open system involving 
effects such as neutrino emission or interaction with a long range electromagnetic field 
whose treatment within a model of the non-relativistic kind studied here is prevented 
by the incompatibility of the necessary Lorentz and Galilean invariance requirements. 

Using the specification (4) in conjunction with the Noether identity (3) we see that 
the dynamical force balance requirement (11) will be expressible simply as 

E/ X M = 0, (12) 

X 

where, for each constituent with label x, the correspond amalgamated force contribution 
is defined by 

h = n+^r-f^,- (13) 
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3. The thermodynamic positivity requirement. 

It is now to be recalled that, according to our preceeding work, each of the or- 
dinary local constituent 4-force densities in the Noether identity (3) will be given as 
the sum of an acceleration contribution and a (gauge dependent) 4-momentum transfer 
contribution by the formula 

n = n + > ( i4 ) 

in which 7r* is the relevant (gauge dependent) 4-momentum covector that is defined in 
terms of the corresponding locally defined material 4-momentum covector and the 
background gravitational field potential by 

^ = /^-ra x <^, (15) 

where m x is the relevant particle rest mass parameter, and the (gauge independent) 
acceleration contribution is specified - using a cross barred symbol - by 

fl = 2»x , (16) 

which therefore satisfies 

<n = o . (it) 

In the particular case of the entropy current labelled by the value x=0, the relevant 
particle rest mass vanishes, and (as illustrated by the example of the historic Landau 
model [3]) the corresponding thermal momentum will be directly identifiable with the 
local temperature covector that is obtained as the partial derivative of the material 
lagrangian density A with respect to the entropy current, 6 M = dA/ds^ 1 (for fixed 
values of the other currents) and from which the temperature in the thermal rest frame 
is obtainable as 

e = -<e M , (is) 

i.e. we shall simply have 

m = O, < = e^. (19) 

This means that, according to (15) the corresponding thermal 4-force density covector 
will be given just by 

f\ = 2s^ [v Q lA + Q tl V v s\ (20) 

with the implication that its time component in the thermal rest frame will be given 
simply by 

uif v = -ev„ s " , (21) 

which will be negative (since the temperature 6 must always be positive) by the second 
law requirement (2) in any system that is closed in the strong sense of being fully self 
contained, so that the external forces f£ t on the right of (11) will vanish. In the more 
general case of an open system, in which there may be heat (i.e. thermal energy) loss 
at a rate (per unit volume in the thermal rest frame) given by 

Q = <flu, (22) 
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due to mechanisms such as neutrino emission, the local formulation of the second law 
of thermodynamics will no longer take simple form (2) but will be given by the more 
ubiquitously valid condition 

Q + 9V,/ > . (23) 

4. Transfusive dissipation 

We use the term transfusion to designate processes by which particles of different 
species are converted into each other by various chemical or nuclear reaction processes, 
which we shall distinguish by capital Greek letters. The elementary process in some 
such reaction, with label S say, will involve the creation of a number N° say of particles 
of the species with label A, using the convention that the early Latin capital A ranges 
over the same values as the late Latin index x except for the zero value reserved for 
the entropy, i.e. it is subject to the restriction A^ 0. (A simple example of such a 
process is the decomposition of a Helium nucleus into a pair of neutrons and a pair of 
protons, so that if we attribute the labels 1,2,3 to neutrons, protons, and Helium nuclei 
respectively, this reaction will be characterised by N~ = 2, N~ = 2, N~ = — 1). In 
any such reaction, the relevant particle creation numbers are restricted to satisfy the 
Newtonian mass conservation condition given by 

J2Nfm A = 0. (24) 

A 

By summing over the rates r s of the relevant reactions, the ensuing particle creation 
rates are obtainable as 

V,<=Er s iV=, (25) 

subject to the proviso that we are dealing with a system that is closed in the weak 
sense [13], meaning that there are no external losses or gains of the relevant particle 
species. Closure in the strong sense, meaning the condition that the model be entirely 
isolated in the sense of being fully self contained, would imply that all the external 
4-force density covectors f* t in (11) should vanish, whereas the weak closure condition 
adopted here corresponds merely to the requirement that the relevant material rest 
frame components should vanish, i.e. 

£X = °- ( 26 ) 

This restriction does not exclude the possibility of the kind of external force that might 
be exerted by interaction with a magnetic field, and since it does not apply to the 
special index value x=0 labelling the entropy, it is also consistent with the possibility 
of a positive value Q > of the heat loss rate (22) due to a mechanism (such as the 
URCA process) of the kind that would necessitate replacement of the simple version 
(2) of the second law of thermodynamics by the more generally applicable version (23). 

In accordance with traditional usage [13] in physical chemistry, it is convenient to 
work with a quantity of the kind for which De Donder introduced the term affinity, 
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according to a specification of the form 

^ = -5^' (27) 

in which the quantities £/S are the relevant energies per particle of the various species 
involved, so that the affinity measure the net energy release in a elementary process of 
the kind (characterised by the label S) under consideration. 

A local static chemical (or nuclear) equilibrium state is one in which the relevant 
affinities vanish. When the deviations from such a state are not too large, it is naturally 
to be presumed that the reaction rates will be linearly dependent on the affinities 
according to a prescription of the form 

r s = E^Al } , (28) 

in which the coefficients n S9 form a symmetric matrix that must be positive in order to 
ensure the positivity of the energy release rate ^2r B A^ e y. This implies that the particle 
creation rates will be given by the formula 

V„< = -£S, fl £f e} , (29) 

B L ' 

in which E AB is a positive indefinite matrix given by 

E AB = ZN-k S9 N*. (30) 

a,* 

This matrix is not positive definite, but only positive indefinite, because it evidently has 
a null eigenvector provided by the set of particle masses, which will satisfy the condition 
Z)2 A m B = by the mass conservation law (24). Other such null eigenvectors will 

B 

be provided by other relevant charge number conservation laws. (The ordinary mass 
conservation law is interpretable as the Newtonian limit of what is given in a relativistic 
theory by the baryon conservation law.) 

The meaning of the foregoing reasonning is unambiguous under conditions of the 
kind most commonly considered in physical chemistry in which there are no significant 
relative motions of the various constituents. However in the more general circumstances 
we wish to deal with here the meaning of the prescription (27) will be affected by the 
choice of reference systems in specification of the energies S^. The most obvious 
possibility is to evaluate the energy in the thermal rest frame, with respect to which 
the energy per particle will be given simply by 

= ~« ■ (31) 

Such a specification does of course depend on the ether frame used for the definition 
of the 4-momentum: it can be seen from the analysis of the preceding article [1] that 
under a transformation characterised by a Galilean boost velocity vector = 7^^V^/9 
it will be subject to a change given by the rule 

£ A i ^ q = SI + m A ^V „/3 - ImH 2 . (32) 
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However the mass conservation condition (24) can be seen to ensure that this gauge 
dependence will cancel out in the corresponding thermal affinity, which is given by 

Al = - E = E N-ny; = Al . (33) 

Although it is independent of the choice of ether frame, the specification (33) of 
the thermal affinity «4| will still be indeterminate in applications for which there is no 
well defined thermal reference frame, a disadvantage that does not apply to what may 
be termed the natural affinity, which is given in terms of the corresponding natural 
energies \\ b Y 

A\ = - E Nlxt = M ■ (34) 

The natural energy is specified for each constituent as the corresponding chemical po- 
tential as measured with respect to its own local rest frame, which means that it will 
be given by 

X$ = = -« + \ mAv l ~ mA( P , (35) 

in which the gauge dependence of the three separate terms on the right cancels out, to 
give 

X A = X A ■ (36) 

In the absence, at this stage, of any clear idea of which, if any of these two (ther- 
mal and natural) alternatives may be most appropriate for general purposes, we shall 
proceed in terms of a compromise using a mixed energy S^ e y that is defined in terms of 
a parameter e by 

£A } = (l- e )£A +eX A > (3?) 

which means that it will transform according to the rule 

£( e} ^ q e} = £( e} + m A (l - e)«V,/3 - \b 2 ) , (38) 

so that as before, in consequence of the mass conservation condition (24), the corre- 
sponding affinity (25) will be invariant, 

A\ e} = A\ e} . (39) 

In this paramatrised weighting scheme the special thermal and natural cases are given 
respectively by 

•^{0} = "4-0 ' -^{1} = • / \ ' (40) 

and the general case we shall have 

Al } = (l-e)AI + eA^. (41) 

When the chemical reaction rates and the relative velocities are sufficiently small 
it will not matter what value is chosen for the parameter e, a consideration that pre- 
sumably accounts for the lack of attention to this issue in the standard literature on 



non-equilibrium thermodynamics [13]. However if the chemical reaction rates or the 
relative velocities are too large, the distinction may become important and in such a 
case the question of what would be most physically realistic would ultimately need to 
be decided on the basis of a microscopic analysis of a kind that does not yet seem to 
have been sufficiently developed. It might turn out that for an accurate out of equilib- 
rium thermodynamic description it would be most appropriate to use something more 
complicated than the kind of weighted mean adopted here. It is to be remarked that an 
affinity of the natural kind, as characterised in our present notation by the weighting 
ansatz e = 1, is what was implicitly used in earlier relativistic work [9, 10], whereas 
use of an affinity of the thermal kind, as characterised by e = 0, was implicit in more 
recent and specialised relativistic work [11]. It will be found below that the latter is 
more satisfactory for applications involving superconductivity. 

5. Viscous dissipation 

In a systematic approach to the construction of phenomenological dissipation laws 
that are consistent with the thermodynamical inequality (23), the first step is to evaluate 
its left hand side by contraction of the 4-force balance equation (12) with the thermal 
rest frame unit vector so as to obtain an identity of the form 

£ Kll , = ? (<^K + <V„r* " - < /*) , (42) 

where is the (purely spacelike) velocity difference between the particular unit flow 
4-vector u£ and the thermal rest frame unit 4-vector u^, i.e. 

<9 = < - < > ( 43 ) 
and is the gauge independent force density contribution given by 

h = h - *f?»< = n + v ^ " - £ m • ( 44 ) 

Let us now generalise this specification to a parameter dependent force density 
contribution given by an expression of the analogous form 

in terms of a combination given for 0, as a function of the weighting parameter e of 
the preceding section, by the formula 

f {e}fi = f A , - {< - e + \rn-v A ^))v v nl , (46) 

so that in particular we shall have 

/((,},= /V (47) 
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We complete the specification for the thermal case with a formula of the rather different 
form 

fly, = fl + E - 6 ( X A M + \m*v A9lt )) V„< , (48) 
in order to obtain a sum over all constituents that is the same as for the ordinary forces: 

E/ { A e} ,= E/*. (49) 

X L ' X 

This specification has been set up in such a way that, unlike the original (canonically 
defined) 4-force vectors which are frame dependent unless the number currents 
are separately conserved, the adjusted 4-force densities given by (46) and (48) will 
have space projected parts that are always unaffected by Galilean (and even Milne) 
transformations, i.e. they will satisfy the invariance conditions 

-rfh* = ^ f U» ■ ( 5 °) 

The motivation for the introduction of the parametrically adjusted force densities 
given by the rather elaborately contrived definition (46) is that it enables us to obtain 
a particularly simple and evocative expression for the entropy term in (42). Using the 
total mass conservation condition that is obtainable from (24) in the form 

E m A V u n u A = , (51) 

A#0 

in conjunction with the restrictions (7) and (26) it can be seen to follow from (42) that 
will be given - for any chosen value of the weighting parameter e - by 

Q + 9VX = -E^ } V,< - ErrV,< - E< ff e}it . (52) 

In view of the second law requirement (23) to the effect that the left hand side of 
(52) should be positive, the choice of admissible dissipation laws will be restricted by 
the condition that it should be such as to ensure that the sum of the terms on the right 
of (52) should also be positive. Although many other (generally more complicted) ways 
involving various kinds of cross coupling are conceivable, it will be adequate for most 
purposes to do this in the most obvious way by ensuring that each of the three sums 
on the right of (52) is separately positive. 

In so far as the first of these terms is concerned, this desideratum of positivity is 
already satisfied by the ansatz of the previous section, which - for the chosen value of 
e - gives an expression of the form 

- E ££ } V„< = E Af e} ^A( e} , (53) 

whose positivity is evidently ensured by the condition that the transfusion matrix k S9 
should be positive definite. 
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To deal with the second term, we exploit the possibility of rewriting the negative of 
the contribution of each constituent as 

rrV,< = 7^ x ^jr. (54) 
where 8£ a is the symmetric spacelike expansion rate tensor given by 

e^ = Y ( ^uK\ e^t a = o. (55) 

This tensor is decomposible in a well defined (Galilean and even Milne gauge indepen- 
dent) manner into tensorially irreducible parts in the form 

0f = < + i* x <r , (56) 

where a£ v is the trace free shear rate tensor and 9 X is the scalar expansion rate, as 
characterised by 

7^ = 0, # x = 7^r- ( 57 ) 

This enables us to write the negative of viscosity term in (42) as 

E rl = E 7mp7-t x "Vf + \ E 7m^ x ^x • (58) 

X X X 

In order to ensure that this total is negative as required, the obvious generalisation 
of the ansatz that is familiar in the case of a single constituent fluid is to postulate that 
the relevant stress contributions are given by an expression of the form 

r x^ = _2E?f Y <" - £C XY # Y 7^, (59) 

Y Y 

where i] XY is a positive definite definite or indefinite but in any case (by the Onsager 
principle) symmetric matrix of shear viscosity coefficients, and ( XY is a similarly sym- 
metric positive definite or indefinite matrix of bulk viscosity coefficients. 

6. Ordinary resistive dissipation 

To complete the determination of the dynamical equations of motion it remains to 
specify the space components of the force on each constituent. The most obvious way of 
doing this in such a way as to ensure consistency with the total force balance condition 
(12) is to take them to consist of sums of pairwise interaction contributions in the form 

7 A "7f e}l , = E/ x ^, (60) 

1 ' Y 

subject to the conditions 

/ XYfi = -/ YxM , V XY,1 = 0. (61) 

It is possible to conceive situations in which a more elaborate construction procedure 
might be needed, but we shall not envisage such complications here. 
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Proceding on the basis of the ansatz (60) we must now consider how the admissible 
forms of the two-constituent interaction force densities / XY M are restricted by the second 
law of thermodynamics. As we have already chosen rules that ensure the positivity of the 
first two terms on the right of (52), this restriction will amount just to the requirement 
of positivity of the final term, which will be given by 

-^<Jm,= Ev.r^^v.r"- (62) 

X 1 ' ^ X,Y X,Y 

in which the (gauge invariant) relative velocity v£ Y and the corresponding (gauge de- 
pendent) covector v XY ^ are defined by 

w £v = <$ ' <0 = <-<> *W = 7^< Y • (63) 

The obvious way to fulfil this requirement is to suppose that the forces are due just 
to resistivity of the ordinary kind, which means that for each pair of distinct constituent 
label values j^Y the corresponding contribution will be given by an ordinary positive 
resistivity coefficient Z XY = Z YX > according to the specification 

/ YX " = ^ Y < Y . (64) 
It follows that the third term on the right of (52) will be given by 

- E < ffo „ = \ E ^ X X Y 7,,< Y , (65) 



which shows that it does indeed satisfy the required positivity condition. 

Having thus obtained an appropriate resistivity formula (60) for the gauge invariant 
force density components 7 Ml fh\ v we can imediately use the defining relation 

h = /{}„ - (n - e W + \m-v A ^))V v nl , (66) 

to provide a corresponding formula for the original unadjusted (gauge dependent) ma- 
terial force density components, which will be given by 

-rh = e ^ ay < a - (yvt - e (y^xt + k mAv A)) v -< • (67) 



6. Superfiuid drag dissipation 

The simple kind of resistive dissipation mechanism described in the previous section 
will not be operational in the case of a constituent that is superfiuid. To deal with 
such cases, let us use indices l,J that range over the values (if any) of x referring to a 
superfiuid constituent, while using indices C,D that range over the values of x referring 
to the remaining - normal - constituents. 

On a sufficiently small - mesoscopic - scale, superfluidity can be dealt with directly 
in terms of the kind of model set up in the preceding section by imposing the following 
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conditions, of which the first and most obvious is simply the requirement that the 
relevant superfluid viscosity and resistivity coefficients should vanish, i.e. 



jyix = o , £™ ; Z IX = , (68) 

which implies, by (59) that that the corresponding viscous tension contributions r lfiu 
will vanish. The next condition is that that the superfluid constituents are not directly 
subject to any external force, 

^ = 0. (69) 
so that we shall be able to make the identifications 

J\ = Jl • (70) 

The final condition that needs to be imposed is that we should be able to make the 
identification 

f{e}n = f{0}fJ. ' (71) 

where (by definition) 

which will follow from (70) subject to the requirement that we adopt the thermal affinity 
ansatz, 

e = 0, (73) 

or else that the relevant chemical reaction rates are set to zero so that the superfluid 
particle creation rates V I/ n 1 ' all vanish, in which case the value of e will not matter. 
When the conditions (68), (69) and (73) are all satisfied it can be seen that the equations 
of motion set up in the preceeding sections will be consistent with the restraint to the 
effect that the relevant superfluid vorticities 

w\ v = 2V[X] , (74) 

should vanish, as necessary for the existence of corresponding local superfluid phase 
scalars y? 1 such that N{K l v = HWuip 1 where N t is the number (2 for the usual case of a 
Cooper type pair) of constituent particle in a boson of the superfluid condensate under 
consideration. 

For application on larger scales such a mesoscopic description is inadequate, and 
must be replaced by a macroscopic description that allows non vanishing superfluid 
vorticities w 1 , which are interpretable a representing the average effect of a fibration by 
microscopic vortex tubes on which the mesoscopic irrotationality condition breaks down. 
This interpretation means that although it does not have to vanish, the macroscopic 
vorticity 2-form of a superfluid constituent can not have an arbitrary algebraic form 
but must satisfy the degeneracy condition 

= » (75) 
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(which means that its matrix rank is not four as in the generic case but two, since 
its antisymmetry excludes the possibility of an odd valued rank) in order for the null 
eigenvectors that generate the vorticity flux two surfaces to exist. As remarked in the 
context of the analogous relativistic problem [11], the only obvious way of setting up a 
force law that satisfies this condition is to postulate that it should have the form 

fh^ri^V", (76) 

in which V" is some vector such that the combination u" + V" is one of the null 
eigenvectors generating the flux surfaces of the vorticity flux zu 1 ^, i.e. such that 

t*W« + V") = o- 

The superfluidity ansatz (76) can be applied within the framework set up in the pre- 
ceding section by supposing that V v is decomposible as a sum of separate contributions 
v from the various non-superfluid contributions, in the form 

V v = Y,V QV , (77) 

c 

so that the corresponding space projected force contributions in the decomposition (60) 
can be taken to be given by 

r^n./V^^. (78) 

As in the preceding section, we now need to find a procedure for choosing the force 
contributions in such a way as to ensure the positivity of the final term on the right of 
(52), which will now be given as a combination of normal and superfluid contributions 
in the form 

- E < /f 0} „ = \ E v CD " + E v cl . (79) 

X C,D C,I 

As before, we can deal with the normal part by taking the relevant force contributions 
to be given for c^D by a set of normal resistivity coefficients Z CD = Z CD > according 
to the simple specification 

/ DCfl = ^%, (80) 

which will automatically take care of the positivity of the first term on the right of (79), 
but the superfluid contributions f ic M will need to be handled in a different manner. 

The remaining positivity requirement that still needs to be satisfied, is that of the 
second term on the right of (79), which can be rewritten as 

E W ICAt = E^^r = E^rw^r > ( 81 ) 

C,I C,I C,I 

where the last step is obtained by substituting (77) in (76) and using the identity 
u if\o}[i = 0- We now procede in a manner analogous to that by which the ordinary 
resistivity forces were introduced above, which means that we ensure the positivity of 
the separate terms in the sum on the right of (81) by taking each vector V cfl to be 
given by an expression of the form 

^^-cO^l-VV^, (82) 
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where c c is a positive drag coefficient and we have included a positive normalisation 
factor given by the magnitude w 1 > of the (gauge independent) space projected 
vorticity vector 

w^ = \e^w\ p , (83) 

according to the specification 

{w^f = i 7 M VMp ro W = l^vt". (84) 

The prescription (84) is evidently equivalent to the adoption of a force law of the 
form 

f l °> 1 = c°n 1 w 1 ±Z u v c , (85) 

where 

Adding up the resulting contributions we finally obtain 

TP m „ = n^w\ p V^ (87) 

which (by the identity u"f l ^ v = 0) automatically provides the required result (76) 
with 

V^-MY^^X- (88) 

To relate this 4-dimensionally covariant formulation to the traditional Newtonian 
terminology using a 3+1 decomposition based on some particular choice of ether frame 
vector e M , it is useful to introduce the (frame dependent) vorticity surface generating 
unit 4-vector u*f and its associated 3-velocity vector vff = — e M by the defining 
conditions 

In terms of such a vorticity flux velocity vector, the degenerate vorticity 2-form c^ 1 
will be expressible as 

^ = (e^ + 2t [fl e u]ap r )w^ , (90) 
and the (rank 2) projection tensor in (85) will be given by 

-L^= it - v?t v - (w')- 2 w^wi . (91) 
It can thus be seen that the projected velocity on the right of (85) will be given by 

_Lim u » = v g _ ( W .)-W't; „- . (92) 



7. The limit cases of convection and pinning 

If some resistivity coefficient, Z XY is very large, the corresponding velocity difference 
will tend to be very small v£ . In such a case it will often be convenient to a use a 
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simplified dynamical treatment based on the relevant convection ansatz to the effect 
that the velocity difference in question should actually vanish. To deal with this kind 
of convective limit, i.e. a case in which the number of independent normal velocities 
is smaller that the number of independent chemical constituents, it is convenient to 
introduce a new kind of index (U) , distinguished by surrounding angle brackets, to label 
the independent velocities and the corresponding comoving subsets of constituents that 
are characterised as equivalence classes by 

ce(u) & «S = « ( V (93) 

In particular we shall use the label (0) for the class of constituents that are convected 
with the entropy, so that 

CG(0) & v^ = 0. (94) 

For each such class of comoving constituents, it will be useful to define combined 
values of additive quantities such as force density and stress using notation illustrated 
by the examples 

fV= E f°, &= E r«2 M = £ r°", (95) 

and more particularly for the ordinary dynamical 4-force density by 

/<?> = E f° = 2u» E n c + E tt° V„n£ . (96) 

ce{u) ^ {U> ce(u) ce(u) ^ 

Using this notation we can regroup the terms on the right of the general entropy 
creation formula (52) in the form 

Q + 6 V„s" = - E %V v n: - E r ^ "V^ - E v£ „ - E vj f\ 0} „ . (97) 
using the notation introduced above, according to which 

t(u) _ f(u) ~{u)v f(u) /qo\ 

The first term on the right of (97) is the rate of chemical energy release which we 
deal with exactly as before in (29) by setting 

V„< = -ES AB £f 0} , (99) 

where H AB is the same reactivity matrix as was introduced above. The second term on 
the right of (97) is the viscous energy dissipation rate which we deal with in the same 
manner as in (59) by setting 

T <u) w = _2 E TyWWfT^ - E ( {u)(v) 0, v)1 ^ , (100) 

(V) {v > (V) { ' 
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the only difference being that the positive shear viscosity coefficients 77 WW = r^ww 
and bulk velocity coefficients ^woo = £00 w now only need to be specified for the 
restricted range of the index (u) labelling the comoving equivalence classes rather than 
for the full range of the normal constituent label C. It will similarly be sufficient to 
specify just a restricted range of positive resistivity coefficients Z< u )< v ) = Z<y)( u ) for 
(u) 7^ (v) to specify the resistivity contribution in a force formula whereby the terms 
in the original expression (62) are regrouped in the form 

Y"fl u L= E/ww + E/^ 1 ", (101) 

1 1 (V) I 

in which the resistive force density terms are given by 

f(u)(v)t* = zw< v >v » . (102) 

J (V)(U) V ) 

It remains to specify the vortex drag terms f l < u ) ^ = —f {u)l , which combine to give 
the space convected superfluid force densities as 

-r/[o}*= s/w, (103) 

1 ' (U) 



so that the last two terms in (97) can be recombined in the form 



According to the reasonning of the preceeding section, these contributions should 
be given by expressions of the form 

for a set of generalised velocity vectors that add up to give a sum 



in terms of which we shall get 



and 



V v = e yw, (106) 

(V) 



^« + ^) = 0, (107) 



The final drag dissipative term in (104) can thereby be rewritten as 

£ tW/./ 1 ™" = S ¥ ( : ) ^,^, (109) 

(u),i y ' (U),l { 1 

By the same reasonning as in the preceding section, we can ensure the required 
positivity of the total (103) by adoption of an ansatz of the form 

v^v = - C iu)( w >)-^ w > v y (u) , (no) 
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which is equivalent to setting 



for a set of positive drag coefficients cj u) where J, like I, ranges over the set of superfluid 
index labels. 

There is however an extreme limit known as pinning - the analogue for a superfluid 
constituent of convection in the normal case - representing what will occur if some 
drag coefficient say is very large, in which case the normal flow vector u £ will 
be constrained to lie in the relevant vorticity surface, a requirement that is evidently 
expressible as the condition 

^;> = °> (H2) 

which is evidently sufficient to ensure that the corresponding dissipation term in the 
sum (109) will simply vanish. 

This can be achieved by taking the term / h <u>a i t Q b e given by the dissipative drag 
prescription of the form (111) that applies to the other contributions / H < v >^ for (v) 7^ (u) 
- and to all the corresponding force contributions for the unpinned constituents - but 
instead by the alternative ansatz 



y ( u)» = v n _ y. VW. (113) 
that is equivalent to the formula 

which is chosen in such a way as to ensure that, according to (103), the total force 
acting on the pinned constituent will be given by 

7^/f 0} „ = ^7^VvV (115) 

For any value of the superfluid constituent index I, not just for the pinned index 
value H we have been considering, knowledge of the contravariant space projection 
7^/joj.y will be sufficient for the complete specification of f\^ v due to the identity 
u if\o}ti wn i cn implies that we shall have 

f\, } v = l^l Vp f\o }p (116) 

with r ) xilv defined as usual by r ) lfJ , v u^ = 0, 7 lM1 ,7^ p = 8£ — t^u^. The formula (115) thus 
implies that for the constituent with label H that is pinned to the set of currents with 
label (u), we shall have 

AV = ™hK„vV - ^>><y > (11 7 ) 
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with the evident implication that u t, ) f{o) ll wm vanish, and hence by the definition (72) 
that the second term in (117) will drop out. Thus we are finally left with the simple 
formula 

ff =n H w* v,\ , (118) 

in which the right hand side is interpretable as a gauge invariant version of the Joukowski 
formula for the Magnus effect. It is easy to see - using the definition (72) and the 
decomposition ufa = + V(u) H ~ that the application of this force law (118) is indeed 
equivalent to the imposition of the pinning condition (112). 

It is to be remarked that the phenomena of pinning and convection are physically 
rather similar in that they both can both b considered as constraints representing 
the effect of extremely strong dissipative coupling. However the way they have been 
dealt with here mathematically is very different. In the case of pinning the dynamical 
equations have been adjusted in such a way that, after having been imposed as as 
an initial value restriction, the relevant restrain will be preserved by the equations of 
motion. On the other hand in the case of convection the constraint has been directly 
imposed at an algebraic level, so as to reduce the number of independent components of 
the system from 4N, where N is the number of constituents (each with its own current 4- 
vector n£) to 4N-3N' where W is the number of independent comotion constraints (each 
of which removes the corresponding 3-velocity components from the list of independent 
variables while leaving the corresponding number density). The most familiar example 
is that of a generic non barotropic fluid, as characterised by just a single independent 
velocity, so that N'=N-1 and the number of independent components is just N+3, 
including as a special case the barotropic fluid model characterised by N=l, N'=0, for 
which the number of independent conmponents redices to 4. 
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Appendix. The basic convective superconducting superfluid model 

In order to set up a large scale neutron star model for the purpose of describing 
the pulsar glitch phenomenon - which is generally recognised to depend on relative 
motion of a superfluid constituent relative to a normal background - the usual kind of 
perfect fluid model will evidently be inadequate. On the other hand it may be hoped 
that a satisfactory description will be attainable without recourse to the very elaborate 
kind of model involving separate allowance for the many degrees of freedom (such as 
those of the electomagnetic field that plays an essential role in the mechanics of the 
external magnetosphere) that would need to be taken into account in a highly accurate 
treatment. As a reasonable compromise, for use in such a context, the following basic 
convective superconducting superfluid kind of model would seem to be appropriate. 

The proposed basic model has three independent constitituents of which one is 
superfluid while the other two are subject to a convection constraint. Thus in the 
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notation of the previous section it is characterised by N=3 and N'=l which means that 
it has 9 dynamically independent components. 

In the context for which it is intended, the three constituent currents are to be 
considered as consisting of an entropy current 

s" = s<, (119) 

a superfluid neutron current 

< = n n <, (120) 

and a normal current 

< = n c <, (121) 
that convects the entropy flux, with which it shares the unit flow 4-vector 

< = <, (122) 

and that is to be interpreted as representing the flux of all other baryons. This normal 
baryon current is to be thought of as consisting not just of protons but at deeper 
levels also of hyperons, whose charge is neutralised by an ambient lepton gas consisting 
not just of electrons but at deeper levels also of muons. In the crust layers the normal 
baryon constituent will also include the neutrons that are confined within atomic nuclei. 
A realistic treatment of the crust (whose outer layers, at densities below 10 11 g/cm 3 , 
contain no unconfined superfluid neutrons at all) would need the use of a model of a more 
elaborate elastic conducting solid kind, whose formulation, in the covariant Newtonian 
approach developped here, will be left for future work. A Newtonian treatment will 
in any case be inadequate for an accurate treatment of the deeper levels, for which a 
relativistic version [11] of the model would of course be needed. 

As in the relativistic version, there will be a baryon conservation law having the 
form 

V,< + V,< = 0. (123) 

In the Newtonian approximation, with which we are concerned here, this law is to be 
interpreted as representing the conservation of rest mass, on the understanding that 
both - superfluid neutron and normal (protonic or other) kinds of baryon are treated 
as having the same rest mass, m say, per particle, while there is of course no rest mass 
associated with the entropy, i.e. we have 

m = O, m c = m, m n = m . (124) 

Within the few per cent level of accuracy - the most that can be expected from a New- 
tonian treatment in this context - this common rest mass m can be chosen indifferently 
to be either the mass of the hydrogen atom or simply the bare proton mass m p (not 
to mention the value traditionally preferred by chemists, which is one sixteenth of the 
mass of an ordinary oxigen atom).. 

In order to characterise a particular model of this type, the essential element - 
which is all that is needed in the conservative case - is the pressure function \1> or 
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equivalently its dynamical conjugate the master function A. nt , whose unconstrained 
version is given as a function of the 4- vectors s M , n£ and in a gauge invariant 
manner - meaning that it can depend only on the the 6 scalars consisting of the number 
densities s n c and n n together with the relative velocity magnitudes v c % v n $ and v nc from 
which the corresponding internal momentum contributions are obtainable as the partial 
derivatives given by the variation law 

6A int = ® v 8s v + X c M: + X a J< . (125) 

The pressure function is obtainable from the master function, and vice versa, by a 
Legendre type tranformation expressible as 

M> = A int - e v s v - x> u c - xlK , (126) 

so that its variation law will have the form 

5V = -s»5& u - <<% - ■ (127) 

The associated stress momentum energy density tensor will be given [2] by a formula 
of the standard form 

T% = s^l + < < + < < + m» . (128) 

in terms of the corresponding set of complete 4-momentum covectors, of which that of 
the entropy is simply 

< = e„. (129) 

while those of the massive constituents are given in terms of the corresponding frame 
dependent 3- velocity covectors v clx — j^u^ and v nfl = ^^u^ by 

^ = X\ + mv c ^- m(\v 2 c + 0)t M , (130) 

^ = X% + mv n ^ -m(\vl + 0)t M , (131) 

where is the gravitational potential. 

In terms of the corresponding thermal, normal, and superfluid vorticity forms, 
namely 

W \v = 2V [M 7r t] » W \v = ^[^l] > W % = WfrK] ■ ( 132 ) 

the associated 4-force covectors will be expressible by the defining formulae 

/J = s^ + e M V^, (133) 
/J = « + W, (134) 

n=<^i,+^<- (135) 

It is to be remarked that if we wanted to describe the superfluid on a mesoscopic 
scale (large compared with the microscopic particle separation lengthscales but small 
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compared with the intervortex spacing) we would need to impose the restraint that the 
supernuid vorticity w^ v should vanish, but that it will in general have a non zero value 
on the macroscopic scale (meaning one that is large compared with the intervortex 
spacing) for which the present treatment is intended. 

If we were dealing with a conducting (as opposed to convective) model the infor- 
mation needed to characterise the dynamical evolution of the system would consist of 
a complete specification of all three of the 4-force covectors that have just been listed, 
which is equivalent to a the specification of the three creation rates V u s u , V u n", V^n^, 
and of the three corresponding space projected 3-force vectors 7^ 'fl 7 A "V£> 7 / "V" (the 
simplest possibility being that of the strictly conservative case for which all three 4-force 
covectors are set to zero). In a convecting model of the kind we wish to consider here, 
we still neeed to specify the creation rates V u s v ', V u n^, V u n^, as well as the space pro- 
jected 3-force vector 7^/° of the superfluid constituent, but in view of the constraint 
(122) we do not need a separate specification for the corresponding thermal and normal 
baryon contributions, but only for their sum 7^ 'ffi where ffi is the combined 4-force 
density defined by 

f { $ = fl + fl- (136) 
The foregoing variational specification of the separate 4-momenta as functions of the 
corresponding currents requires that the master function be defined not just for con- 
vectively constrained configurations but even when there is a relative motion between 
the entropy current and the normal baryon constituent. Such a general, unconstrained 
specification will indeed be available if the convective (9 component) model under con- 
sideration has been obtained as a high resistivity approximation from an unconstrained 
(12 component) conducting in which all three constituents move independently. It will 
however be more economical from a mathematical point of view to avoid the introduc- 
tion of redundant information from such an unconstrained ancestor model, and to work 
entirely within the framework of a reduced variational framework in which the master 
function A. nt and the associated pressure function \1/ are specified only for the range of 
variables allowed by the convectivity constraint (122). 

In such a reduced formulation, the 9 independent components can be taken to consist 
of the 4 components of the superfluid current vector n^, the 4 components of the normal 
baryon current vector n", together with just the density s of the entropy current, whose 
velocity is not independent but given by that of the normal baryon current. The most 
general variation that is allowed within this reduced formulation will be given by an 
expression of the form 

M int = -6 5s + x f 5n» c + X * 5< , (137) 

which does not provide a specification of the separate thermal and normal baryon 
momentum covectors @ M and Xu, but only of the scalar temperature G and the amal- 
gamated normal momentum covector x^% that can be evaluated within the ancestral 
unconstrained framework as 

e = - u ;e u , x { $ = x c + -(et^ + ej. (138) 
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The corresponding complete amalgamated normal momentum covector 

^ = ^ + -f(e^ + e,)> (139) 

will be obtainable in the reduced formulation as 

7r<g = X <J + mu^ - m(^ c 2 + 0)^ . (140) 

It is to be emphasised that it is possible for physically different ancestor models - 
as characterised by different values of M and x M f° r given values of the independent 
currents - to engender the same reduced model, in the sense of providing the same 
values for and an d hence also for the pressure function (127) and the stress 
energy tensor (140), which will be expressible as 

= A int + Qs - X { i < ~ XIK , (141) 

and 

T1 = < 7r< c j - QsH v + < nl + . (142) 

The formalism of the reduced formulation is not quite so elegant, but it has the avantage 
of avoiding the introduction of operationally redundant information singling out some 
particular one of the compatible unconstrained ancestor models. In the framework of 
the reduced formulation the combined force (136) will be expressible as 

/<£> = 2n c "V[„7r<g + 7r<g V„n c " + sV,& - ^V,(0^) . (143) 

We thereby obtain an expression of the form 

rV { jg„ = 2 7 ^<V [i ,7r< c > p] + s 7 ^V M + (f"(x { H - X\) + m<)V,< , (144) 

for the corresponding adjusted (gauge invariant) 3-force density vector, whose analogue, 
for the free neutron current, will be given simply by 

TTo}, = 7 W « P , < P = 2V[,7r n p] . (145) 

In order to complete the determination of the dynamical equations of the model, it 
is necessary to choose the rules specifying the values of these space projected 3-force 
vectors 1^ f^\ u and ^ v f^ u , and of the two independent creation rates V v s u and V v n^ 
- of which the latter, by (123), determines V v n". The simplest possibility is of course 
that of a strictly conservative model for which the forces and creation rates all vanish. 
What we want to consider here is the more general case in which there is internal 
dissipation by the mechanisms described in the preceeding sections and perhaps also a 
non vanishing heat loss rate Q (due to neutrino emission) but we shall suppose that 
although the system may thus be open in the thermodynamic sense it is nevertheless 
isolated in the sense there are no external contributions to the space projected force 
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densities, nor to the time component of the force acting on the free neutron current, so 
that in the notation of the preceding section we shall have 

&v = °> & = Qt v . (146) 

According to the general principles developed in the preceding sections, the space pro- 
jected forces will therefore be given by expressions of the standard form 

-TfUv = /< c > n "-V v r< c >^, (147) 

r u f{o } v = -f {c)nf "-V l/ T^». (148) 

In a high temperature version of the model, the mutual interaction force density f^ n ^ 
would be given in terms of an ordinary positive resistivity coefficient Z^ n by an ex- 
pression of the form 

j(c>n M = Z (c)n v £ . ( 149 ) 

What we are particularly interested in here however is the low temperature version of 
the model, in which the free neutron current is a superfluid, which means that instead of 
being given by an ordinary resistivity formula of the form (149) the mutual interaction 
will be given by a vortex drag formula of the kind given by the formula (111). 

In the present case, the vortex drag force density will be given by the expression 

/ (c)v = -n»f</i c)p , (150) 

in which is the positive drag coefficient, and the vector V^f 1 ^ is defined by 



K <c> " = -^r-T^x > = h^ pa <P<« > ( 151 ) 



W 1 

or alternatively just by 

^ C> " = «<S„> ( 152 ) 

in the non-dissipative large limit case of vortex pinning. In any such (pinned or 
unpinned) superfluid model, the free neutron current will not be subject to any viscosity 
force, i.e. we shall have 

r n ^ = 0, (153) 

but for the combined thermal and normal baryon current contribution there will in 
general be a non vanishing viscosity contribution of the standard form 

r (c> [iv = _2yf) a ^ _ £<c>0 c7 v» ? (154) 

in which rf^ and are positive shear and bulk viscosities (of which the latter will be 
negligibly small in many circumstances) and and 6 C are the trace free and trace 
parts of the normal constituent's expansion tensor as given, according to (56) by 

9? u = 7 CT(/1 V CT <) = a? + \e c ^ v . (155) 
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To complete the specification of the model it remains to give the prescription for 
the creation rates. In the present case the only relevant kind of generalised chemical 
reaction is one whereby a free neutron is created by a processe such as "dripping" out 
of a confined state within a neutron star crust nucleus or by inverse beta decay of a 
proton in the core, so that the corresponding creation numbers will be N n = 1 and 
N c = —1, and the corresponding chemical affinity will be given by the formula 

-4{0} = £{0} ~~ £{W ' ( 156 ) 

in terms of the relevant thermal rest frame energies as defined in terms of the relevant 
(ancestral or reduced) 4-momentum covectors by 

% = = ~^ {C J , % = ■ (157) 

The superfluid particle creation rate can thus be seen to be given in terms of the 
relevant transfusion coefficient k by an expression of the standard form 

V u n^ = KA m , (158) 

in which the affinity is given in terms of the relative flow velocity magnitude v nc by the 
manifestly frame independent formula 

-4{0} = x\ ~ X$ - \rnvl c , (159) 

where 

x\ = ~<xl = -<xl c> , xl = -< Xl = ~<xl ■ (160) 

The last thing we need to complete the specification of the model is the value of the 
energy emission rate Q that determines the entropy rate via the formula (97), which 
gives 

Q + ev^ = kA\ $} + ^ liivlpa + C <c> ^c + ^K <c> "7m,K <c> " ■ (161) 

Cri 

The complete set of dynamical equations for the 9 independent components (those of 
the space vectors and v£ together with the scalars n c , n n and s) is thus completed: 
it consists of the creation formulae (123), (158) and (161), together with the pair of 
3-force equations (147), (148) (as made explicit by the prescriptions (150), (151, (153), 
and (154) for drag and viscosity). 

It is to be remarked that final term in (161) will drop out in the large df> limit 
for which the drag prescription (151) is replaced by the vortex pinning prescription 
(151). In a similar way the first term in (161) in the large k limit for which the system 
will be maintained in a state of chemical equilibrium as characterised by the condition 
-4{0} = which will have the effect of reducing the number of dynamically independent 
components from 9 to 8. 
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The particular case of a thermodynamically closed model - which may be a good 
approximation for processes occurring on a short timescale - will be obtained by setting 
Q = 0. However it will often be more realistic to take Q to have a value that is positive 
and monotonically increasing as a function of the temperature 9 to allow for losses by 
URCA type neutrino emission processes. For processes ocurring over a sufficiently long 
timescale the temperature sensitivity of Q near some emission threshold value may be 
sufficient to justify the use of a simplifying approximation whereby the temperature 
9 is held fixed at the threshold value in question, thereby determining the value of 
s and hence of the creation rate V^s^. This will reduce the number of dynamically 
independent components from 9 to 8 (or, in the chemical equilibrium case, from 8 
to 7) so that (161) will no longer be needed as a dynamical equation of the system, 
but will merely serve for the purpose of calculating the corresponding value of Q in 
case it might be needed. A simple extreme special case of such a fixed temperature 
thermodynamically open variant of the model is the zero temperature limit for which 
9 and s both vanish. 
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